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P> ■ Abstract 



We present a derivation of the effect of the classical field configuration to the diffusion 
equations. Using the formalism of the thermo field dynamics we propose a systematic 



Oh! 

^ ■ and consistent way to treat the classical background and to calculate it's effect at any 



order of the perturbation theory. We treat the classical field in genuine quantum field 
theory formalism as a coherent state and show how the propagators and self energies 



H are altered m this case. 
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In statistical physics, the non-equihbrium phenomena are traditionally described 
by the Boltzmann equations. Classically speaking, these equations describe the trans- 
portation of particles in the background of the statistical fluid, disturbed by the col- 
lisions between particles. The Boltzmann equations can be derived in a variety of 
ways, the most fundamental being from the Dyson-Schwinger (DS) equations of the 
relativistic, statistical field theory. 

If the non-equilibrium state is caused by spontaneous symmetry breaking (SSB), a 
spatially and temporally non-homogeneous classical field describes the propagation of 
the phase transition. In that case we would like to have a way to calculate the effects of 
this SSB field to equilibration of the particle gas. This can be done by working out the 
derivation from the equations of the quantum field theory to the Boltzmann equations 
and keeping in mind that the symmetry breaking field, i.e. the order parameter, is 
a macroscopic classical field. From the Boltzmann equations one could derive e.g. 
diffusion-like equations for the particle densities. However, in the present paper, to 
make contact with the studies of the Refs. H, ^, we by-pass the derivation of the 
Boltzmann equations and write the diffusion equations directly from the DS equations. 

The proper way to deal with classical fields in quantum field theory is to consider 
them as coherent states This means that the field quanta propagate coherently, 
in constructive interference, and form a macroscopic field. This state conserves its 
coherence and the expectation value of the field evolves according to classical equations 
of motion as far as the number of quanta in it remains large 0. In this proper 
treatment all phenomena, e.g. decay of the coherent state, are naturally incorporated 
to the theory. 

Recently there has been debate, how a classical field configuration, describing the 
propagation of the symmetry breaking phase transition, should be treated and what 
is its effect, especially in the context of the particle gas out of equilibrium in the 
early universe ||, It seems not to be clear, how to incorporate the properties 
of the domain wall into the quantum field theory. Various methods from more or 
less heuristic |l| to a kind of semi-classical ones 0, has been adopted leading to 
somewhat different results. 

A well-argued technique to account the effects of the classical field configuration 

is of especial importance when considering the baryogenesis during the electroweak 

phase transition [^. Near the boundary of the broken and the symmetric phases the 

particle gas is in non-equilibrium state, thus fulfilling one of the necessary Sakharov 

conditions for the baryogenesisQ. The phase boundary acts then as a separator of 

quantum numbers. By various methods the quantum numbers are converted to baryon 

^The other conditions, baryon number violation and CP- violation, can also be fulfilled: CP- 
violation through the usual CKM-matrix based mechanism (or other mechanism in the case of the 
extensions of the Standard Model ) and baryon number violation due to the sphaleron processes 
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number transported to the broken phase and recognized as the baryon asymmetry of 
the universe. 

We aim to show in the present letter how the classical field, which describes the 
symmetry breaking phases in different regions in space, can be adequately accounted 
using coherent states and thermo field dynamics (TFD) [|l^. The coherent state is the 
counterpart of classical field in the quantum field theory. On the other hand, TFD 
gives the way to treat the non-equilibrium thermodynamics consistently |]1T . 



In the formalism of TFD, to account the doubling of the degrees of freedom of 
the finite temperature field theory, one introduces the operators S,k-, ^k, ix^ 
where K = {ko, k) labels the energy and the momentum. These operators obey the 
commutation relations 

= {2ixY6^'\P - K), = {2ixY6^'\P - K) (1) 

and vanishing other commutators. These are related to the operators a^^, ax, cl^k 
and ajf, which are used to construct the Heisenberg operator fields and their tilde 
conjugate^, by the Bogoliubov transformation 




B 



ax 

T 



-ax 



B-\ (2) 



Here the Bogoliubov matrix B has a unit determinant and leaves the form ^^'^ — 
1^^^ invariant. Thus it is not uniquely defined. Intuitively, a simple choice is the 
so called symmetric representation or a = 1/2 -gauge, when is the hermitian 
conjugate of (which is not true in general). This means that the Bogoliubov 
matrix satisfies the relation B^^ = t^B^t^, where = diag(l, —1). The benefit of the 
symmetric representation is that the left- and right-handed thermal vacua (see below) 
are conjugates. However, in the symmetric gauge there appears to be some problems 



with the fermionic case [12 



Another common choice for B is the so called a = 1 -gauge |]T3[, which we shall 
adopt in the present paper. Then 

bh''_V";i. (3) 

where n is the bosonic number density. The benefit of this gauge is that also fermions 

can be treated consistently. Note further, that in finite temperature the energy and 

^ The finite temperature Hamiltonian is unbounded from below[p2[. Tlie ordinary fields represent 
the particle-hke excitations with positive energy while the tilde conjugated fields represent hole-like 
excitations with negative energy. 
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momentum labels ko and k are not related by a dispersion relation. This is a conse- 
quence of the broken Lorentz-invariance in the thermal bath. 

The Fock-space associated with the ^-operators has the right-handed vacuum state 
\R) , which is annihilated by the operators and ^k] and the left-handed vacuum 
state (L|, which is annihilated by the operators and ^f^. The one-dimensional 
realization of the vacua, which can be readily adapted to the more general case, is (in 
the a = 1 -gauge) 

\R) = exp [fa^d^) |0) 

{L\ = (0|exp(a5), (4) 

where |0) is the vacuum state of a and d -operators and / = e~^'^ is the statistical 
weight. With these definitions, statistical averages of observables correspond to the 
expectation values, 

TT{pd[A,A^) = {L\d[a,amR) , (5) 

where p = exp{—/3ujN) is the statistical operator and N = A'^A is the number operator. 

At this point we make a generalization to the formalism of TFD by shifting the 
operator A with a parameter t] so that the new statistical operator is 

= exp -r]){A- r])) , (6) 

where bar denotes complex conjugate. The corresponding right-handed vacuum of 
TFD is 

\R,7]) =exp(/a^5"f + (l-/)(r7a"^ + r/a"^-|r/H) |0) (7) 

and left-handed vacuum is left unchanged: {L,ri\ = {L\. Note that the shift in 
generalizes to SSB field theory where one shifts the symmetry breaking field with the 
vacuum expectation value. In the zero temperature limit the state (0) reduces to 

\r]) = exp (rja'^ + f/a^ — jr/p^ |0) . (8) 

Hence the expression (|^) can be regarded as a finite temperature generalization of a 
coherent state This is in accordance with treating the symmetry breaking field as 
a coherent state. 

The ^ operators, which annihilate the vacua \R,f]) and (^^,''71 are now defined 
through 





(9) 
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With these definitions, it is easy to calculate that (we are still considering the 
one-dimensional case) 

. {L,r]\a^a\R,ri) ^ 

where 

i.e. number density is a sum of the contributions of the background "field" rj and 
quanta generated by the operator A. 

Changing to the field theory, one similarly postulates the left- and right-handed 
vacua, which are annihilated by the raising (^|^ and ^|^) and lowering (^x and fx) 
operators, respectively. These are related to the a-operators by the generalized Bo- 
goliubov transformation (|^) for every mode separately: ax — t]k ~^ C,k- Here rjx is the 
Fourier transform of the background field (f){x), to be discussed below. 

The Heisenberg fields are constructed in two steps. First, one introduces the so 
called generalized free field 4>k{x) fl^ 0, 

M^) = k-'l ^5(|p|-A:)(e-*'*+^P-ap + e^'=«*-^P-al,)U=.„ (12) 

and corresponding tilde conjugated field cpx- These are hermitian fields, but generaliza- 
tion to more complicated cases is straightforward. In the second step one decomposes 



the interacting Heisenberg field A(x) as [|T4| 



SO that the field Ak has a limit lim(_>_oo Ak = (pK in Lehmann-Symanzik-Zimmermann 
(LSZ) sense. Note that the positive definite function Z is a generalization of the wave 
function normalization factor of the vacuum field theory. In principle it is fully deter- 
mined by the canonical commutation relations and the LSZ-condition However, 
in practise its calculation is a tremendous task. 

Perturbation theory is derived now from the finite temperature form of Gell-Man- 
Low -formula (GML) [T^: as usual, the Feynman rules are derivable from it. In the 
presence of the coherent state rix it generalizes for the two point function to 

iD{x, x') = (L, 7]\ T[A{x)A^{x')] \R, r]) = (L, r/| T[^x)^\x')U] \R, r]) , (14) 
where 

U = Texp( r dtHiit)) (15) 



■^Our notation for is slightly misleading here, because it depends on energy and square of the 
momentum k. 
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with the interaction Hamiltonian Hi = H[^]—Hq[^]. Here the thermal Hamiltonian is 
expressed as H[^] = H — H, where H is the ordinary hamiltonian governing the time 
evolution of the Heisenberg fields and H is its tilde conjugate. For tilde conjugation 



rules, see Ref. pO |. 



Thermal doublet notation was introduced in Eq. (|14[): 



H^) = I ~u \ 1 > (16) 



XI 



the field having the same form as the interacting field A, but replaced by (px- 
The interacting thermal doublet field A is defined with a formula similar to (|I6D. 

The zeroth order propagator (which is 2 x 2 -matrix) can be calculated even without 
specifying the interaction Hamiltonian. By a straightforward application of Eqs. (^, 
(11), (HI) and (ll) one obtains 

D(°)(x, y) = Dix, y\=o - i^c{y)Ux)J = D^'\x, x\=o + 5D{x, x'), (17) 



where (/)c is the classical field 



[X 



and J is a constant matrix 



The -D^=o part of the propagator (p!7D has the same form as the usual TFD zeroth order 
propagator without coherent state. However, it depends on r] through the functional 
dependence of Z{K) on rj. 

Note, that in order to use DS equations (to derive quantum transport equations) it 
is necessary to assume that the zeroth order propagator is on mass-shell, in the sense 
that 

(□, + m2)D(°) (x, x') = S^^^ {x ~x')l, (20) 

where 1 is 2 x 2 unit matrix. This implies that rjK and Z{K) are proportional to the 
mass shell delta function S{K'^ —rn?). In what follows, we adopt this simple, although 
generally illegitimate, approximation. This approximation is physically well motivated 
if the particle mass is large enough compared to its (thermal) width, i.e. F <^ m. This 
makes it reasonable to speak about particles because the free path 1/F is long enough 
so that they can be regarded as "free" between interactions. Note also, that because 
particle current can be written in the form 

m^)) = = -\ \{d^-d^)D'\x,y)\_, (21) 



x=y 



*Here we consider complex scalar field with corresponding generalizations in the notations. 
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the particle number N{x) in zeroth order reads 

Nix) ^ (Jj(x)) = / -0^,nM + ^0e(x)U0c(x), (22) 

where mass shell form of Z(K) was used, so that cuk = [k"^ + m^)^/^. Thus the Eq. 
( p^ ) produces the conventional result. 

Within the above framework, it is possible to derive the quantum transport equa- 
tions in the presence of the SSB field. Note that the dynamics of the interacting fields, 
i.e. the Heisenberg equations, are not affected by inclusion of the thermal coherent 
state. The thermal state (whether the coherent state included or not) contributes to 
the propagators only through the boundary conditions for the propagator. 

Although the quantum transport equations (i.e. Boltzmann equations) can be 
derived from the DS equations^ using conventional procedures, we choose to step 
directly to the diffusion equation, which is in many cases sufficient. To begin with, 
the DS equation for the scalar propagator matrix reads 

D{x,x') = D^^\x,x') + Jd^yJ d'^zD{x,y)U{y,z)D^^\y,x'), (23) 
where 11 is the self energy (which can be calculated perturbatively) 



/ nil , 
n = I n- ' " ^^^^ 



There is also a similar equation to Eq. (^) with D and D^'^^ interchanged in the 
right hand side. Note that for the simplest case, constant background field (pc, the 
resumming of mass (in the 0^ -theory) follows easily from Eq. (|23|). 

Expressing the particle current as in Eq. (|21|) and using the DS equations we 
obtain for the four divergence of the current 

^0(a:) + V-j<^(x) = -Jd^yeix'~y'')[u'\x,y)D'\y,x)-Il''ix,y)D^\y,x) 

- D'\x,y)U'\y,x) + D'\x,y)U'\y,x)], (25) 

where {J^{x)) = {n^{x) , j^{x)) . In deriving Eq. ( [25| ) the conditions D^^ + D"^"^ = 

L)12 + 1)21 ^11 ^ ^22 ^ _]-[12 _ ]-[21 ^g^g ^gg^_ 

Next we write corresponding equation for the fermion field ip. The fermion propa- 
gator is defined through 

S{x,x') = -t{T[^{x)^{x')]), (26) 



^For example, see Ref. |13[|. 
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where \1/ is the thermal doublet 

^ix) = y ^= (^(x) - zj'^ix)) . (27) 

The left- and right-handed vacuum states are generalized in an obvious way to be 
the vacua with respect to the fermionic raising and lowering ^-operators. Similarly as 
in the scalar case, using DS equations, one obtains for the divergence of the particle 
current (J^), 

( J^(x)) = inxfrnx)') = Tr [rS{x, yri=y, (28) 

where trace is taken over Dirac indices, 

n4^) + ^-i4^) = -Tr J d'yOix' ~y')[j:'\x,y)S'\y,x) -J:'\x,y)S'\y,x) 

- S'\x,y)J:'\y,x) + S'\x,y)E'\y,x)], (29) 

where S is the self energy of the fermion field \E'. Similarly as in the scalar case, we 
used the conditions 5" + S^^ = S^^ + and S^^ + = -S^^ - S^i. 

Next step would be to express the self energies in terms of the propagators (using 
perturbation theory) and obtain self-consistent equations for the particle currents. We 
do not perform this step in all its length, but choose to study only the effect of the 
inclusion of the coherent state to the continuity equation. 

Because we do not consider fermion fields to have any coherent part, the only direct 
contribution to Eq. (|29| ) appears through the scalar propagator 6D{x,x'). We write 
the self energies in the form 

n = n^=o + and S = S^=o + (30) 

where all coherent state dependence, i.e. dependence on SD, is included in 511 and 
SH. These quantities can be calculated perturbatively from the Feynman rules. The 
bosonic continuity equation reads 

n0(x)+V-j^(x) + Jd^y[u'\x,y)D'\y,x)-U''{x,y)D''{y,x) 

+ D'\x,y)U'\y,x)-D'\x,y)U'\y,x)] ^ 

= 7<A,c«(a;), (31) 

where the coherent field contribution reads 

7<^,cz(x) = -ld'y[6U'\x,y)D'\y,x\=o + 'nlloix,y)SD'\y, 



X 



6n''ix,y)Df^,iy,x) - nli,ix,y)6D'\y, 



x] 



+ 6D'\x, y)Ull,iy, x) + D]^^,{x, y)5Yi^\y, x) 

- 5D'\x,y)nf^,{y,x)-D':^^,{x,y)5n^\y,x)]. (32) 
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The corresponding formula for the fermionic equation is 

n^(x) + V-j^(x) + Ti Jd^y[E'\x,y)S'\y,x)-J:^\x,y)S'\y,x) 
- S'\x,y)J:'\y,x) + S'\x,y)J:'\y,x) 

where 



r)=0 



(33) 



Tr I d^y [5T?\x,y)S^\y,x\=, + T?^^o{^,y)5S^\y,x) 



(5S {x,y)S ^^{y.x) 



21 



x,y)SS {y,x 



6S''{x, y)T}^Uy^ x) - S%{x, y)6T}\y, x) 
5S'\x,y)J:f^,{y,x) + S^Ux,y)SJ:'\y, x) ] 



(34) 



The left hand sides of Eqs. (^Tj) and (^) are the usual terms of the diffusion equation 
after applying Pick's law j = —DVn. In the linear approximation, the last term in l.h.s 
produces the damping term — F n. The right hand sides of these equations implement 
the effect of the background field (pc- 

One can apply Feynman rules to obtain perturbative expressions to the integral 
terms in Eqs. (|3l| ) and (|33D . For example, lowest non- vanishing perturbation calcu- 
lation apphed to ( p^ ) reproduces exactly the result of Ref. 0. In these papers CP- 
violating effects are considered at electroweak domain wall. There the classical field 
is inserted to the theory as an interaction with fermion field through a Yukawa cou- 
pling gYip(f)ip thus producing CP-violating current at two insertion level. Indeed, the 
same result can be obtained by noting that at one-loop level with Yukawa-coupling 
the correction to the fermion self energy reads (a, (3 = 1, 2) 



SJ:{x,yr^ = -gl.S'{x,yr^6D{x,yr^, 



(35) 



and 6S = 0. The right hand side of Eq. 



reads in the one loop approximation 



7^^) = -9y J d'yO{x' - y') [M^Hliy) - MvHli^ 



xTr 



(36) 



This leads to the formula as given in Ref. 0. 

In the present letter we have proposed a method for consistent treatment of clas- 
sical fields in the framework of thermo field dynamics. In the method, a classical field 
is accounted as giving a coherent contribution to the thermal vacuum of the TFD. By 
using perturbation theory, the effects of the classical field can be calculated system- 
atically. In the lowest non-trivial order the calculated contribution of the background 
field to the fermionic diffusion equation reproduces the result of Ref. p|. Because the 



8 



inclusion of the coherent state to the thermal vacuum does not affect the formalism 
of the quantum field theory, all relevant aspects of the underlying model are auto- 
matically incorporated. In particular, possible CP, baryon or lepton number violating 
effects of the model are naturally accounted. This is important e.g. in cosmological 
considerations, where classical field configuration, i.e. coherent state, represents the 
regions of broken symmetry. We will return to these aspects in a forthcoming paper. 
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